Tunnelling dynamics of Bose-Einstein condensate in a four wells loop shaped system 
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We investigated the tunnelling dynamics of a zero-temperature Bose-Einstein condensate (BEC) 
in configuration of four potential wells arranged in a loop. We found three interesting dynamic 
regimes: (a) flows of matter with small amplitude, (b) steady flow and (c) forced flow of matter for 
large amplitudes. The regime of quantum self- confinement has been studied and a new variant of it 
has been found for this system. 



PACS numbers: 03.75.Kk, 03.75.Lm, 74.50.+r 



I. INTRODUCTION 

The behavior of Bose-Einstein condensate in two wells 
potential based on two state approximation has been in- 
vestigated in a number of theoretical and experimental 
papers 0, i, i, i, i, 0, 0). Here we expand this 
method to a four-well system with periodic boundary 
conditions. The main aim of this analysis is to study 
possible ways in which to achieve mass transport around 
a loop and persistent currents. 



II. THE GROSS-PITAEVSKII EQUATION 
WITH THE FEYNMAN ANSAZ 

The behavior of a Bose-Einstein condensate at low 
temperature is accurately described by a nonlinear 
Schrodinger equation, known as the Gross-Piteavskii 
equation (GPE), obtained from the two bodies interac- 
tion Hamiltonian by neglecting the quantum fluctuations 
of the bosonic field. 

The GPE equation that describes a BEC trapped in a 
potential V{r) is: 



~dt 



2m 



+ V(r) +^ I ^ 



(1) 



with a coupling constant g — AnTi^a/m^ where a is the 
scattering length of atom (that has mass m), taken here 
to be positive. 

We take the normalized solution of the time- 
independent GPE for the ith non interacting well to be 
and look for an approximate solution of the form 
(ll|,|9|,ll0|): 



<i{r,t) = Y,^i{t)^i{r). 



(2) 



Substituting this form into Eq.([T]) and writing the 
time-dependent function ipi as ^/Nle^^^ we obtain, after 



an integration over the spatial variables, a system of four 
coupled complex equations: 



ihdt^2 = {El 
ihdt^s = {El 
ihdt^^ = {El 



(3) 



with the condition Xli=i = ^t- Where A^r is the 
total number of particles, Ef is the ground state energy 
for the ith well, Ui is the self-interaction energy of the ith 
well and Ki is a parameter that characterizes the overlap 
between wells. 

It is useful to recast Eq.(|3l) as a system of eight real 
equations: 



hdtNi 



+Ki-i^Ni-iNism{0i - 9i-i] 



+Ki^Ni+i/NiCos{ei+i 

+ Ki_iy^Ni_i/NiCOs{ei - 



(4) 



-UNi-Ef 



with z = 1 to 4 and all the arithmetic on the index 
modulo 4 {N3 = Ni etc.). 

As pointed out in [1] the total number of atoms Nt 
is constant but, in order to have a coherent phase de- 
scription, the phase fluctuations must be small, giv- 
ing a lower bound on the number of particles, namely 
Nt > Nmin — 10^ [11]- One may notice that typically 
the value of /7i is largely independent from z; we thus drop 
the index i and consider all the Ui to be the same. In 
the rest of the paper we will consider only the case of all 
the Ki modulated around the same mean value and the 
case of all the Ki fixed to a certain positive value or set 
to 0. We will call K this mean or fixed value and define 
the two quantity uor and A, such that K = fiuoR/2 and 
A = UNt/2K. The first parameter has the dimension 
of a frequency and we will take its inverse as time unit 
while the second is a dimensionless quantity regulating 
the behavior of the system. 
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FIG. 1: The lay out of our system with the couphng constants 
evidenciated. 
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FIG. 2: (Color online) The lines A,B,C and D are the plots 
of population versus time in the linear regime for the first, 
second, third and fourth well respectively. In (a) the mass 
surplus increases during the resonant driving from time t = 
to t = T, where r = A/ujr. For t > r there is a dissipa- 
tionless flow around the loop. In (b) the modulation is not 
interrupted and the amplitude beats are visible. In these sim- 
ulations Nt = 10^ and A = 500. 



III. SMALL AMPLITUDE REGIME 

Starting with an arbitrary small population inbalance 
in one of the four wells we can in principle, by mod- 
ulating the coupling constants, amplify it and make it 
spin around the ring. In the case of symmetric wells 
{Ef = Ej for all i, j) we simulated the system with cou- 
pling constants of the form: Ki = K{1 -\- {—ly sm{wt + 
0)), i.e. a periodic oscillation at the frequency w = 

V^UNtK ^2k'^, that is the resonance frequency of 
the Eq.dlj) linearized around the value Ni = Nt/^ and 
Oi — 6i-^i = for all i. The result of our simulation is 
shown in Fig. [21 starting with an arbitrary small popula- 
tion excess in one well the resonant driving increases the 
population difference and makes it spin around the loop. 
We are effectively causing the system to spin by modula- 
tion of the coupling constant in the same way one could 
make a ball turn in a dish by raising and lowering its 
edges at the right moment. The direction of the spinning 
is imposed by the initial phase (j) of the perturbation rel- 
ative to the position of the initial imbalance. 

If we stop the modulation at t = r we can see 
in Fig. [2] (a) that the particle surplus keeps going 
around the loop and will only be dumped by dissipa- 
tion ( 0, [12I, ll^, llSl ) . However the dumping time is 
of the order of some substantially longer than the 

time-scale we are interested in. We can thus neglect it. 
In any case we must notice that this process can prove 
to be very difficult to measure in practice, due to the ex- 



tremely small amplitude of the oscillations. Indeed one 
can not amplify the imbalance over a certain small per- 
centage of the average well population as in this case, due 
to the change of resonance frequency of the system with 
the amplitude of the oscillations, we see periodic beats 
of the amplitude ( Fig. [3](b)). As A decreases, the peak 
amplitude of these beats increases but their frequency 
decreases. This means that, in order to have beats large 
enough to be observed, their period will be of the order 
of some tens of I/oor and thus the dissipation could not 
be a negligible problem anymore, due to the intrinsically 
longer time-scale implied. As stated in [12] the applica- 
tion of a periodic modulation of the potential well can 
stabilize the tunneling dynamics against dissipation and 
thus the beats in amplitude could be effectively observ- 
able. In any case, the exact study of these dissipation 
effects is out of the scope of our present work. 



IV. NONLINEAR REGIME 
A. Phase imprinting 

In the nonlinear regime, exploiting the fact that the 
phase is defined modulo 27r, we can have a constant de- 
phasing between neighboring wells. If we set all the pop- 
ulations at the same value ^A^t, all the Ef and all the Ki 
at the same constant values and all the phase differences 
between adjacent wells to 7r/2, building up in this way 
a phase difference of 27r around the loop, we have that 
in Eq.(|4|) the gain and the loss term for each well equi- 
librates each other, and the relative phases of the wells 
stay constant. Being in this case the system invariant 
under rotations of 7r/2 and being the evolution of the 
populations dependent on the relative phases, we can ex- 
pect that the only possible behaviors of the system can 
be a steady flow of particles around the loop or no flux at 
all. In both cases, from our simulation we will see noth- 
ing but the fact that all the populations stay constant at 
the value of \Nt. 

In order to probe what is happening we simulate an 
abrupt cut of the loop at one point (by setting one cou- 
pling constant to zero for t > r) and observe at the re- 
sults, shown in Fig. [3l If effectively there was a steady 
flow we would expect it to continue by inertia so that 
the last well before the cut would be filled with atoms 
while the others would be depleted; this is exactly what 
emerges from Fig. [3] (a). In experiments it is often quite 
difficult to control the exact value of the coupling con- 
stants and also to assure that all of them have the same 
value. A coupling constant different from the others will 
behave as a bottleneck, thus creating an oscillation of 
populations superimposed on the steady flux (Fig. [3] 
(b)). In order to perform any experiment it would then 
be necessary to assure that these oscillations do not hide 
the current we are interested in. In any case this turns 
out not to be a major problem as both the period and 
the amplitude of these oscillations increase with A. It is 



3 




Self-confined well 

\/v\/w^^A/\/\yw 



(b) 



30 ^^^^^T^CXTv^^CTv^^ 



Time (1/co ) 



Time {1/co ) 



FIG. 3: (Color online) (a) Populations versus time plot for 
a system with a phase difference of 27r around the loop. For 
t < 0.5/oor there is a steady flux and the populations stay 
constant, for t > 0.5 /ur the link between two wells is cut 
and by inertia the condensate tends to flow toward the last 
well before the cut, finally for t > 2/oor the inertial effect 
is exhausted and the condensate starts to flow back. The 
lines A,B,C and D are relatives to the first, second, third and 
fourth well respectively, (b) System with a persistent current 
in presence of a coupling constant different from the others. 
The population of one well is plotted for several values of this 
coupling constant (Ka = K^Kb — 1.2K,Kc = IAK^Kd = 
1.6K). In these simulations Nt = 10^ and A = 100. 



then a priori possible, knowing that we have a given de- 
gree of incertitude in the coupling constants, to tune the 
value of UNt in order to have enough time to perform 
our measurements before large oscillations develop. 



B. Self-confinement 

As in the two- well case ([l?]) the nonlinearity of the 
GPE leads to the possibility of self-confined states. 

As a condition for a self-confined state to arise we can 
take that in the time in which the phase difference be- 
tween two wells varies between and tt their population 
difference doesn't change sign. Using Eq.(|4]) we can de- 
rive, by means of very rough approximations, an approx- 
imate condition for this to happens: 



2 ,-3^3, ^ ^ 
— tan +1>0. 
3n ^ 4An ^ 



(5) 



Where n is the normalized population imbalance be- 
tween the well in which we are considering the self- 
confinement and the population of an adjacent well, sup- 
posed to be We can check that this condition 
is satisfied both for positive and negative values of n. 
The first case corresponds to the usual quantum self- 
confinement, the second one corresponds to a state of 
self-depletion, in which one well almost empty at the be- 
ginning remains almost empty. Numerical simulations 
show that condition stated in Eq.([5|) works effectively 
quite well. For A = 100 and Nt = 10^ it gives us a 
minimal number of atoms in one well for being in the 
self-confined regime of 31750 while the simulated value 
is 35000 and a maximal number of atoms in the consid- 
ered well for being in the self-depleted regime of 18250 
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FIG. 4: (Color online) Populations versus time plot of (a) the 
self- confinement and (b) self-depletion regimes for A = 100 
and Nt = 10^. Due to the symmetry of our system two 
populations have the same value and so only three lines are 
visible. 



where the simulated value is 15000. It is important to re- 
mark that in the second case the populations of the other 
three wells are smaller than the lower limit for the quan- 
tum self-confinement, so effectively we are considering a 
new effect, not simply a self-confinement of the three full 
wells. We can have an insight into this phenomenon if 
we interpret the usual quantum self- confinement of the 
two-well system as a self-confinement for one well and a 
self-depletion for the others. In the two- well system the 
conditions for having them coincide (obviously if one well 
is empty the other has to be full) but in the four wells 
case they don not coincide anymore and we can thus ex- 
perimentally distinguish them (Fig. [3|). In both cases 
we will have to face the standard dissipation phenomena 
we encounter in the usual quantum self-confined regime 
( [I2I , [Tsj ) but also in this case the time-scale we are in- 
terested in to observe the phenomenon are substantially 
smaller when compared with the typical dissipation time. 



C. Full condensate spinning 

Exploiting the self-confining nonlinearity we can make 
the entire condensate move all the way around the loop. 
Let us suppose that the condensate is at the beginning at 
one well. Increasing the coupling constant toward an ad- 
jacent well, enough in order to make the well exit the 
self-confined regime, causes a flux toward this second 
well. By lowering the coupling constant before the fiux 
changes sign, we find ourselves with the second well in 
the self-confined regime. Repeating this process we can 
make the condensate move around the loop (Fig. [5]) in a 
sort of quantum version of the conveyor belt mechanism 
([16]). The number of atoms in the depleted wells has 
to be quite small (of the order of few percents of Nt) 
in order to avoid having significant spurious population's 
oscillations that could bring the system off resonance af- 
ter very few turns. The process has revealed itself to be 
almost independent from the initial phases of the four 
condensates. This spinning of the whole BEG is possible 
because we have no intensity resonance, as we had before, 
and so the resonance frequency remains the same. In any 
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FIG. 5: (Color online) The condensate is moved around the 
loop by making the different wells to switch between the nor- 
mal and the self-confined regime. The lines A,B,C and D are 
the populations respectively of the first, the second the third 
and the fourth well. After two full turn the modulation is 
stopped and the spin ceases. In this simulation Nt = 10^ 
and A ^ 100. 



case, because of the long period implicated (of the order 
of some centains of 1/ujr), dissipation could play an im- 
portant role. However, contrary to what happens in high 



frequency oscillation in one self-confined well, here we 
have only very low frequency population's fluctuations. 
Thus the only significant dissipation source should be the 
spontaneous atom loss that is normally possible to keep 
relatively small for the time-scale we are interested in. It 
is not straightforward to calculate the correct modulation 
frequency analytically but simple numerical integration 
of Eq. (|3|) is sufficient to find it accurately enough for any 
experiment. We may notice however that contrary to 
what happens in the others two cases, as soon as the 
modulation is interrupted, the spinning stops. 

D. Conclusion 

We have shown that, even though the dynamics of BEC 
in two-well system is a well studied problem, the exten- 
sion to four wells allows us to highlight new and interest- 
ing aspects. It also gives us the possibility to understand 
better well known phenomena. We would like to stress 
that we used the four wells configuration simply because 
it is the less complex nontrivial loop configuration, but 
that the phenomena we highlighted have nothing specific 
to the number four, being all of them trivially expand- 
able to loops made up of a larger number of wells. S.D.L. 
would like to thank Yvan Castin for useful comments and 
corrections. 
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